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Исследована возможность применения алгоритмов муравьиных колоний для реализации криптоанализа 
шифров перестановок, результатом применения которых к открытому тексту является криптограмма, полу- 
чаемая путем перестановки символов открытого текста в определенном порядке. Показано, как эта пробле- 
ма может быть сведена к классической задаче о назначениях, решаемой с помощью алгоритма муравьиных 
колоний. Приведен алгоритм решения, дан пример работы муравьиного алгоритма. 
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Введение. В последние годы интенсивно разрабатывается новое научное направление с назва- 
нием «природные вычисления», которое объединяет математические методы, содержащие прин- 
ципы природных механизмов принятия решений [1]. Как отмечено в [2], научное направление 
«природные вычисления» объединяет такие разделы, как эволюционное программирование, 
нейросетевые вычисления, алгоритмы роевого интеллекта, муравьиные алгоритмы, генетические 
алгоритмы. В [3-6] рассматривались методы организации криптографических атак на традици- 
онные симметричные криптосистемы, использующие шифры перестановки и замены, а также на 
блочные криптосистемы с использованием методов эволюционной оптимизации и генетического 
поиска. 

Однако, как отмечено в [7], структуры генетических алгоритмов являются «слепыми» по- 
исковыми структурами с присущим им рядом недостатков. Поэтому представляет интерес приме- 
нение конструктивных эвристических методов, идеи которых заимствованы у живой природы или 
физических процессов и в которых решение задачи строится поэтапно путем добавления нового 
компонента к частично построенному решению. К методам данного вида относят и муравьиные 
алгоритмы, основу которых составляет имитация самоорганизации муравьиной колонии. Как от- 
мечено в [1], в противоположность примитивному поведению отдельных агентов поведение всей 
колонии оказывается достаточно разумным. 

Отметим, что муравьиные алгоритмы исследуются с середины 90-х годов, и на сегодняш- 

ний день известны их применения к задаче о коммивояжере [8], квадратичной задаче о назначе- 
ниях [9], задаче о раскраске графа [10, 11], задаче маршрутизации в коммутационных сетях [12], 
задаче маршрутизации транспортных средств [13]. В данной работе мы рассмотрим возможный 
подход для реализации криптоанализа шифров перестановки и покажем, как эта проблема может 
быть сведена к классической задаче о назначениях, решаемой с помощью алгоритма муравьиных 
колоний. 
Понятие шифров перестановок. В качестве первичного признака, по которому производится 
классификация шифров, используется тип преобразования, осуществляемого с открытым текстом 
при шифровании. Если буквы открытого текста при шифровании только меняются местами друг с 
другом, то данный шифр относится к классу шифров перестановок. Отметим, что основные виды 
шифров перестановок описаны, например, в [14, 15]. В общем случае результатом применения 
данного класса шифров к открытому тексту является строка символов (криптограмма), получае- 
мая путем перестановки символов открытого текста в определенном порядке. 
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Таким образом, полученная криптограмма включает только те символы, которые состав- 
ляют открытый текст. Отсюда следует, что задача определения открытого текста заключается в 
определении позиций для назначения символов криптограммы таким образом, при котором целе- 
вая функция, определяющая оптимальность исходного текста, достигает экстремума. То есть дан- 
ная задача криптоанализа, по сути, является частным случаем задачи о назначениях, цель кото- 
рой — определить экстремум затрат, необходимых для обмена ресурсами между всеми объектами. 

В соответствии с [16] задачу о назначениях сформулируем в следующем виде. Определим 
Х,Е1, если объект 1 назначен в пункт /, и Х,=0 в противном случае, С; — затраты на передачу 
объема ресурсов из пункта г в пункт /. В этом случае оптимизационная модель запишется сле- 
дующим образом: 


п п 


В=>. > С‚Ху — экстр. 


НЕ 


при ограничениях 


> ХХ, #=1,2,...,П, 
+ 


и 
ре 


где и — число объектов и мест их размещения. 

Применительно к задаче криптоанализа будем полагать, что С, — вероятность того, что за 
символом в позиции 7 должен следовать символ в позиции 1+1, кроме этого введем параметр О, 
показывающий, насколько фрагмент текста из г символов носит осмысленный характер, т.е. сов- 
падает с словарным запасом языка. В этом случае оптимизационная модель будет иметь вид 


п п 
К р ОС,Х; — шах. 

Отметим, что элементы С, задаются в виде матрицы размерности ихи (п — число симво- 
лов текста). 

Таким образом, множество вариантов решений определяется числом перестановок Р=и! 

без повторений и символов, входящих в шифртекст в и позициях. Комбинаторный характер этой 
задачи приводит к необходимости использования метаэвристических алгоритмов. 
Алгоритм решения. Таким образом, общее значение целевой функции К, получаемое в каждом 
конкретном варианте назначения символов в позиции, может быть аналогично [9] определено как 
длина маршрута, соединяющего выбранные элементы декартова произведения [номер позиции, 
номер символа], т.е. как 


Очевидно, маршруту с большим значением А должна соответствовать более высокая концентра- 
ция феромона РГ, которая используется в качестве вероятности выбора очередного маршрута, 
представляющего очередной вариант назначения символов в позиции, новыми муравьями- 
агентами. 

Отметим, что в соответствии с [1] любой муравьиный алгоритм независимо от модифика- 


ций должен быть представим в следующем виде: 

1) создание популяции муравьев; 

2) поиск решения; 

3) обновление феромона. 

Отметим, что на этапе 1 выбор стартовых точек для размещения популяции муравьев за- 
висит от ограничений, накладываемых условиями задачи, так как для каждой задачи способ раз- 
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мещения является определяющим. Либо муравьи размещаются в одной точке, либо в разных с 
повторениями, либо в разных без повторений. На этом же этапе задается начальный уровень фе- 
ромона, который представляет собой небольшое положительное число, чтобы вероятности пере- 
хода в следующую вершину не были нулевыми. Определение вероятности перехода между сосед- 
ними вершинами и нового уровня феромона после обновления может производиться по форму- 
лам, приведенным в [1]. 

Таким образом, в соответствии с [1] и [9] алгоритм включает следующие этапы. 

1. Случайным равновероятным образом выбираются т вариантов маршрутов, и вычис- 
ляются значения целевых функций Ку, К›,...Ви. 

2. Комбинациям #1 размещения символов К в позиции г присваивается весовой коэф- 
фициент 

ЛеЕКЬ [1,2,....т В (1) 


3. Для каждой комбинации {# вычисляется результирующая концентрация 
Ва-УЛы . (2) 


Для тех комбинаций (4, которые ни разу не встретились в выборке т, задается нижнее значение 
концентрации феромона 


Еши=а лтахркл (3) 
где 0<а<1. 
4. В соответствии с формулой, приведенной в [17] — 
(д=т(о--р), (4) 
производится имитация испарения феромона со всех комбинаций {, по которым прошли 
муравьи. 


5. После определения нового количества феромона производится возврат муравьев в 
начальные позиции и определение вероятностей размещения символа А в позиции Г в новом 


маршруте 
Ра=Ещ/ Е) (5) 


Совокупность указанных вероятностей образует матрицу вероятностей размещения ихя. 

6. В соответствии с вычисленными вероятностями Р‚. формируется 4т новых маршру- 
тов (4<1), для которых определяются критерии К.и,...К„-„, И далее производится выборка из 
т лучших вариантов. Если оптимальное значение критерия не изменяется в течение достаточно 
большого количества циклов, то поиск завершается с найденным значением К„„, в противном 
случае длина пути обнуляется и производится возврат к шагу 2 алгоритма. 

Отметим, что, так как шифртекст может содержать повторяющиеся символы, то будем по- 

лагать, что циклы в маршруте не запрещены и длина маршрута ограничена числом симво- 
лов шифртекста. 
Демонстрационный пример. Рассмотрим функционирование представленного выше алгоритма 
на демонстрационном примере. Пусть задана строка символов Б К СОА. Требуется определить 
возможную перестановку символов, входящую в словарный состав языка. Вначале на основе сло- 
варя русского языка составим матрицу С.;, показывающую вероятность того, что за символом 1 
может следовать символ / (рис.1). 
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Б К С [@) А 
0,01 0,01 0,1 0,5 0,6 


0,01 0,01 0,01 0,5 0,4 
0,05 0,08 0,05 0,6 0,3 
0,6 0,3 0,5 0,02 0,1 
0,6 0,6 0,6 0,1 0,01 


рообяп 


Рис.1. Матрица С, элемент С; которой определяет вероятность 
соседства в тексте символов { и] 


Определим количество муравьев т=5 и поставим их в соответствие каждому символу. 

Итерация 1. 

1. На 1-ом этапе выберем случайным образом т маршрутов, представляющих варианты 
размещения символов в позиции, и определим значения их критериев К;, ..., А„. Пусть выбраны 
следующие маршруты: 

















1. КСБАО К=0,01+0,05+0,6+0,1=0,76 
2. БКОСА №&=0,01+0,5-0,5+0,3=1,31 
3. СКАБО Вз=0,08+0,4+0,6+0,5=1,58 
4. ОАКБС К4=0,1+0,6+0,01+0,1=0,81 
5. АКБОС &5=0,6+0,01-+0,5+0,5=1,61 














Поскольку 1, 2, 4, 5 варианты далеки от словарного состава языка, то умножим их на ве- 
совой коэффициент О=0,5, а 3 вариант, наиболее близкий к словарному составу, умножим на ко- 
эффициент О=0,9. Получим следующие значения критериев: А,=0,35; Ю=0,655; Юз=1,422; 
К.=0,405; К5=0,805. 

2. На 2-м этапе всем комбинациям размещения символов в позиции присваивается весо- 
вой коэффициент в соответствии с формулами (1) и (2). Для тех комбинаций, которые не встре- 
тились в выборке из т маршрутов, зададим нижнее граничное значение концентрации 
Ети=а тпах/и.,. Значение а определим как а=0,1, тогда Р.„=0,14. Матрица результирующих кон- 
центраций феромона будет иметь вид, показанный на рис.2. 


Позиции 
| 2 3 4 5 
Б 0,655 0,14 1,155 1825 0,14 
5: 0,35 2,88 0,405 0,14 0,14 
3 С 1,42 0,35 0,14 0,655 1,21 
о 0,405 0,14 0,655 0,805 1/72 
А 0,805 0,405 1,42 0,35 0,655 


Рис.2. Матрица результирующих концентраций феромона после 1-й итерации 


3. Далее в соответствии с формулой (4) проведем испарение феромона, определив в со- 
ответствии с [17] р=0,6 (рис.3). 
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Позиции 
1 2 3 4 ь 
0,262 0,14 0,462 0,73 0,14 
0,14 1,152 0,162 0,14 0,14 
0,568 0,14 0,14 0,262 0,484 
0,162 0,14 0,262 0,322 07 
0,322 0,162 0,568 0,14 0,262 


Символы 
рооятп 


Рис.3. Матрица концентраций феромона после его испарения после 1-й итерации 


4. Вычислим вероятности размещения символов в соответствующие позиции по формуле 
(5) (рис.4). 


Позиции 
1 2 3 4 о 


0,18 0,08 0,29 0,46 — 0,08 
0,09 0,66 0,1 0,09 0,08 
0,39 0,08 0,09 0,16 0,28 
0,11 0,08 0,16 0,2 0,4 

0,22 0,1 0,36 0,09 0,15 


Символы 


рообяп 


Рис.4. Матрица вероятностей размещения символов в позиции после 1-й итерации 


5. В соответствии с вычисленными вероятностями сформируем новые маршруты, вы- 
брав &/=0,8. Разместим муравьев в следующие позиции символов в соответствии с вероятностями 
размещения: С, А, Б, О. Случайным образом сформируем 4 маршрута: 




















1. АОКБС К‹=0,1+0,3+0,01+0,1=0,51 
2. СКБОА К-=0,08+0,01+0,5+0,1=0,69 
3. БКАОС К;=0,01-+0,4-0,1+0,5=1,01 
4. ОКАБС Ко=0,3+0,4+0,6+0,1=1,4 


Для вариантов 6, 7, 8, 9 определим значения О’ соответственно: 0,5; 0,9; 0,6; 0,5. 

Получим: 
К=0,255; В=0,621; Юз=0,606; Ю.=0,7. 

Из полученной популяции выберем 5 вариантов с лучшими значениями целевой функции. 
Это будут варианты: А>, Аз, Кь, Ю7, №. Обозначим: А:=А›=0,655; Е›=Аз=1,422; Вз=Аь=0,805; 
К.=КВ7=0,621; В5=Ко=0,7. 

Итерация 2. 

1. Матрица результирующих концентраций феромона будет иметь вид, показанный на 
рис.5. 

2. После испарения феромона получим матрицу концентраций, показанную на рис.6. 

3. Матрица вероятностей размещения символов в позиции показана на рис.7. 
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Позиции 
1 2 3 4 > 


0,655 0,14 1,426 2,122 0,14 
0,14 4,203 0,14 0,14 0,14 
2,043 0,14 0,14 0,655 1,505 
0,7 0,14 0,655 1,426 1,422 
0,805 0,14 2,122 0,14 1,276 


Символы 
рообятпт 


Рис.5. Матрица результирующих концентраций феромона после 2-й итерации 


Позиции 
| 2 3 4 5 
Б 0,262 0,14 0,57 0,848 0,14 
2 к 0,14 168 014 0,14 0,14 
8 с 0,82 0,14 0,14 0,262 0,602 
9 0,28 0,14 0,262 0,57 0,56 
А — 0;322 0,14 0,848 014 0,51 


Рис.6. Матрица концентраций феромона после его испарения после 2-й итерации 
Позиции 

1 2 3 4 5 

0,15 0,06 0,29 0,45 0,07 

0,07 0,76 0,07 0,06 0,07 

0,4 0,06 0,07 0,13 0,3 

0,15 0,06 0,13 0,3 0,29 

0,18 0,06 0,44 0,06 0,26 


Символы 
рооятпт 


Рис.7. Матрица вероятностей размещения символов в позиции после 2-й итерации 


4. Случайным образом сформируем 4 маршрута, разместив муравьев в следующие 
позиции: С,С,А,Б. 

















1. СКБОА К‹=0,08+0,01+0,5+0,1=0,69 
2. СКОБА ®=0,08+0,5+0,6+0,6=1,78 
3. БКАОА &=0,01+0,4-0,1+0,1=0,61 
4. АКББА К№-=0,6-+0,01-+0,01-0,6=1,22 








Для вариантов 6, 7, 8, 9 определим значения О соответственно: 0,9; 1; 0,5; 0,5. Получим: 
Кв=0,621; К›=1,78; Кз=0,305; Ко=0,61. 
Далее вновь из полученной популяции выберем 5 лучших маршрутов с лучшими значе- 
ниями целевой функции. Это будут варианты: Ау, К›, Юз, К, К. Обозначим: Ю:= В7=1,78; 
К›=ЕК=1,422; Вз=ЁВз=0,805; Ю.=Кь=0,7; К-=Ав=0,621. 
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Итерация 3. 
1. Матрица результирующих концентраций феромона будет иметь вид, показанный на 
рис.8. 
Позиции 
1 2 3 4 5 


0,14 0,14 1,426 3,902 0,14 
0,14 5,328 0,14 0,14 0,14 
3,823 0,14 0,14 0,14 1,505 
0,7 0,14 1,78 1,426 1,422 
0,805 0,14 2,122 0,14 2,401 


Символы 


роохлапй 


Рис.8. Матрица результирующих концентраций феромона после 3-й итерации 
2. После испарения феромона получим матрицу концентраций, показанную на рис.9. 


Позиции 
1 2 3 4 5 
0,14 0,14 0,57 1,56 0,14 
0,14 2,13 0,14 0,14 0,14 
1,52 0,14 0,14 0,14 0,602 
0,28 0,14 0,712 0,57 0,56 
0,322 0,14 0,84 0,14 0,96 


Символы 
рБроохлхаоая 


Рис. 9. Матрица концентраций феромона после его испарения после 3-й итерации 
3. Матрица вероятностей размещения символов в позиции показана на рис.10. 


Позиции 
1 2 3 4 5 
0,05 0,05 0,23 0,63 0,06 
0,05 0,80 0,06 0,05 0,06 
0,66 0,05 0,06 0,05 0,25 
0,11 0,05 0,30 0,22 0,23 
0,13 0,05 0,35 0,05 0,40 


Символы 
Боохая 


Рис. 10. Матрица вероятностей размещения символов в позиции после 3-й итерации 


4. Продолжая процесс далее, определим 4 случайных маршрута, размещая муравьев в по- 
зиции С, С, С, А. 























1. СКАБО К‹-=0,08+0,4-0,6+0,5=1,58 
2. СКОБА ®=0,08+0,5+0,6+0,6=1,78 
3. СКАБА &:=0,08+0,4+0,6+0,6=1,68 
4. АКОБА Ко=0,6+0,5+0,6+0,6=2,3 


Для вариантов 6, 7, 8, 9 определим значения О соответственно: 0,9; 1; 0,9;0,9. Получим: 
Кв=1,422; К7=1,78; Кз=1,512; К.=2,07. 
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Выберем 5 маршрутов с лучшими значениями целевой функции. Это будут варианты: К», 
К7, Вы, Кз, В. Обозначим: Ю!=К=2,07; Е›=К7=1,78; Юз=Ю!=1,78; Ю4=Аз=1,512; А5=А№=1,422. 

Итерация 4. 

1. Матрица результирующих концентраций феромона показана на рис.11. 


Позиции 
1 2 3 4 5 
0,14 0,14 0,14 8,564 0,14 
0,14 8,564 0,14 0,14 0,14 
6,494 0,14 0,14 0,14 0,14 
0,14 0,14 5,63 0,14 1,422 
2,07 0,14 2,934 0,14 7,142 


Символы 


роодхлп 


Рис.11. Матрица результирующих концентраций феромона после 4-й итерации 


2. Матрица вероятностей размещения символов в позиции после испарения показана на 
рис.12. 


Позиции 
1 2 3 4 5 
0,03 0,03 0,03 0,88 0,03 
0,03 0,88 0,03 0,03 0,03 
0,69 0,03 0,03 0,03 0,03 
0,03 0,03 0,60 0,03 0,15 
0,22 0,03 0,31 0,03 0,76 


роохап 


Рис.12. Матрица вероятностей размещения символов в позиции после 4-й итерации 


3.Разместим 4-х муравьев в позициях С, С, С, А и определим 4 случайных маршрута: 


1. СКОБА № =0,08+0,5+0,6+0,6=1,78 
2. СКОБА К7=0,08+0,5+0,6+0,6=1,78 
3. СКОБО Аз=0,08+0,5+0,6+0,5=1,68 
4. АКОБО Ко=0,6+0,5+0,6+0,5=2,2 


Для вариантов 6, 7, 8, 9 определим значение 0 как 1; 1; 0,9; 0,8. Получим: 
К=1,78; К7=1,78; Кз=1,512; Ко=1,76. 

Выберем 5 маршрутов с лучшими значениями целевой функции. Это будут варианты: 1, 
>, Вз, Кв, В. 
Выводы. Была рассмотрена возможность применения алгоритма муравьиной колонии для реше- 
ния задачи криптоанализа шифров перестановок. Комбинирование алгоритма муравьиных коло- 
ний с элементами генетических алгоритмов, а также распараллеливание алгоритма муравьиных 
колоний может существенно повысить эффективность использования данных методов и вероят- 
ность нахождения оптимального решения. 
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\У.А. РАТКНТ, А.5. ЗЕКСЕУЕУ 


АРРЫСАТТОМ ОЕ АМТ СОГОМУ АЕСОВТТНМ РОК ВЕАМ?АТТОМ 
ОЕ ТКАМ$РО$ТТТОМ СТРНЕВ$ СКУРТ АМАЕУ$Т$ 


Ро5 Ру оЁ аррИсайоп оГапЁЕ соопу а/дог те Гог геайгайоп оЁ {те {тапзро$оп с|рВег$ сгурапа/Уз/5 15 5иед. 
[65 аррИсайоп Ю {те р/ат {ехЕ гезий5 т {те сгурюдгат гесемеа Бу $ Ипа зутро/5 ог {те рат {ехЕ т а сейат 
огаег. [+ 15 пои Вои/ #5 ргоМет сап Бе ге тсеа Юю а са5$/са/ ргоШет о! а!осайоп 5оме4 Бу апЁ соопу а!до- 
пёрт. Тре зоивоп адопЁт 15 мет. Ап ехатрЕ оЁ ап а/догЁйт регогтапсе 15 дезспреа, 

Кеу игога&: сгурапа/з65, рго Мет о! а!осайоп, апЁ а/догЁйт, рреготопе, гапзроШоп с/рВег. 
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